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ABSTRACT 

We show that the usual picture of supersoft X-ray binary evolution as driven by 
conservative thermal-timescale mass transfer cannot explain the short orbital periods 
of RX J0537.7-7034 (3.5 hr) and IE 0035.4-7230 (4.1 hr). Non-conservative evolu- 
tion may produce such periods, but requires very significant mass loss, and is highly 
constrained. 
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1 INTRODUCTION 

The importance of episodes of mass transfer in a semide- 
tached binary on a thermal timescale has recently been em- 
phasized in a number of contexts. Such episodes can arise in 
either of two ways: 

(i) the donor star fills its Roche lobe while already un- 
dergoing thermal expansion across the Hertzsprung gap, or 

(ii) the donor fills its Roche lobe with a mass ratio 
q = M2/M1 (M2 denotes the donor mass, Mi the accretor 
mass) large enough that the Roche lobe radius 7?l shrinks on 
mass transfer more rapidly than the thermal equilibrium ra- 
dius Rtc{M2) (where this exists) appropriate to the donor's 
current mass Af2. 

In case (i) the thermal-timescale episode will end once 
the donor attains a new thermal equilibrium radius, e.g. 
at some point after reaching the Hayashi line. In case (ii) 
the donor is continually trying to expand thermally beyond 
_Rl in order to reach Rtc, so mass is transferred on a ther- 
mal timescale while this condition holds. Generally this case 
involves shrinkage of the orbit to some minimum Roche 
lobe size, followed by orbital expansion as the mass ratio 
q reverses. The thermal-timescale episode ends only when 
Rtc < Rl (cf. Fig. pi; see also van den Heuvel 1992 and 
references therein). Thereafter normal mass transfer contin- 
ues driven either by systemic angular momentum losses (de- 
screeising Ri^) or nuclear evolution of the donor (increasing 

Rtej- 

Case (i) arises in the formation of intermediate-mass X- 
ray binaries with black-hole accretors, such as GRO J1655- 
40 (Kolb et al., 1997; Kolb, 1998). (In similar systems 
with neutron-star accretors, such as Cyg X-2, the condi- 
tion for case (ii) may hold simultaneously - King & Rit- 
ter, 1999; Podsiadlowski & Rappaport, 2000, Kolb et al., 
2000; Tauris et al., 2000). In this paper we are mainly 



concerned with case (ii). This has received most attention 
in connection with the supersoft X-ray binaries (van den 
Heuvel et al., 1992). Thermal-timescale mass transfer from 
a donor initially on or close to the main sequence, on to a 
white dwarf accretor, offers a way of driving accretion rates 
Ml ~ 1O~^M0 yr~^ high enough to allow steady nuclear 
burning. As well as potentially explaining the observed su- 
persoft systems, this process allows the white dwarf mass 
Ml to grow. If one can arrange that Mi reaches the Chan- 
drasekhar mass Mc — 1.44M0 this suggests a way of mak- 
ing Type la supernovae. However the difficulty of computing 
mass transfer on these timescales meant that early studies of 
this process simply used the assumption that mass transfer 
occurred on a thermal timescale, and were therefore unable 
to predict the evolution of the binary parameters (masses, 
period, mass transfer rate). Detailed calculations of these 
have only recently begun to emerge (Deutschmann, 1998) 
but are not yet exhaustive. 

Discussion of the evolution of supersoft X-ray bina- 
ries would be greatly eased if observation provided reliable 
masses. However this is very difficult for several reasons. 
For example Greiner et al. (2000) suggest rather low masses 
Ml ~ O.6M0,M2 ^ O.35M0 for the short-period system 
RX JQ537.7-7034 we shall discuss extensively in this paper. 
But these and similar estimates use the assumption that the 
donor is close to the main sequence. By definition this cannot 
be true in thermal-timescale mass transfer, since the star is 
not in thermal equilibrium; although it may be quite close to 
its main-sequence radius, it can also be considerably smaller 
than this (cf. Deutschmann, 1998). Moreover Greiner et al.'s 
mass function uses the Hell emission line, and thus may not 
reflect the dynamical motion of the accretor. 

Since mass information is so hard to come by, a cru- 
cial test for the thermal-timescale model for the supersoft 
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Figure 1. Schematic diagram with i?tc (full line) and /Jl (broken 
line) versus M2 for Case A mass transfer, after Giannone, Kohl 
&; Weigert (1968). While the Roche lobe curve varies significantly 
with initial WD mass Mn and the specification of mass transfer 
{rj, j3), Rte depends only very little on the evolutionary state along 
the main sequence (indicated by the central hydrogen content 
Xc). Thermal— timescale mass transfer may occur until the two 
curves intersect which happens at different M2 depending on a 
particular choice of Rl- 



X-ray binaries is provided by the discovery of systems with 
fairly short orbital periods. In particular RX J0537. 7-7034 
(Greiner et al., 2000) has a period of about 3.5 hr, and 
IE 0035.4-7230 (= SMC 13) (Schmidtke et al., 1996) has a 
period of 4.126 hr. Since the initial mass ratio Qi must be 
^ 1 for this type of evolution (see below), one expects ini- 
tial donor star masses M2i ^ IM© for a typical white dwarf 
accretor, and thus an initial orbital period Pi ^ 10 hr (see 
Section 3 and Fig. bl) . Evidently considerable orbital shrink- 
age would be needed for such systems to reach the periods 
of RX J0537.7-7034 and IE 0035.4-7230. This is unlikely if 
the mass transfer is conservative, as we shall show. How- 
ever it is probable that much of the transferred mass is not 
accreted by the white dwarf, but blown away from it as a 
wind (e.g. Li & van den Heuvel, 1997), allowing greater or- 
bital shrinkage. Mass loss induced in some way by the mass 
transfer process is probably the only way of significantly in- 
creasing the orbital shrinkage, as other angular momentum 
loss processes such as magnetic braking or gravitational ra- 
diation generally take place on timescales far longer than 
the 10^ yr characteristic of mass transfer in the supersoft 
X-ray binaries. 

In this paper we consider thermal-timescale mass trans- 
fer and investigate how much mass loss from the accretor is 
required if RX J0537.7-7034 and IE 0035.4-7230 are prod- 
ucts of the standard picture of the supersoft binaries. Our 
method is to compute the minimum orbital period analyti- 
cally for specified rates of mass and angular momentum loss 
from the binary. 



2 ORBITAL EVOLUTION 

We consider the orbital evolution of a semidetached binary 
in which a fixed fraction 1 — r; of the mass transferred from 
the donor (star 2) is lost from the accretor (star 1) with 
/3 times the specific angular momentum of the latter. (The 
quantity 1 — r; is called a by King & Kolb, 1995 in their 
general treatment of such 'consequential angular momentum 
loss' [CAML] mechanisms; the use of -q allows more compact 
formulae in what follows.) Following the general method of 
van Teeseling & King (1998), setting their quantities Mu,2 = 
Jsys = 0, /?i = /3, we find 

M^ui = {l-V)M2, (1) 

so that 

Ml = -M2 + M^i = -■nM2- (2) 

These definitions give 

with M = Ml + M2 . Kepler's 3rd law links angular momen- 
tum J and period P of the orbit as 



r = G 



,MfMi P 



M 2-k' 
The evolution of the orbital period then obeys 

P_ 3M2 , 3M2 2(1 - r?)M2 , ,,, „,,^ ,, M2M2 



(4) 



(5) 



M2 Ml M ' """^ ' MMi 

From (|) we have Mi = C - r?M2, M = C+{l-ri)M2, with 
C — Mii -\- riM2i, where Mii,M2i are the initial values of 
the two masses. Thus (|5|) integrates to 

P oc M2^M;^^^-''^-'"'^Ud^-^^ (6) 

for 77 7^ 0, 1, the corresponding expression for r; = being 
given by the limit of this expression as 77 -^ 0, while we get 
the well-known result 



p oc M2"^M^^ 



(7) 



for conservative mass transfer, i.e. t) = 1. In terms of the 
mass ratio q = M2/M1 we can express the relation between 
Af2 and Mi as 



gA'/i ^ M2^ + -{A4li 
V 



leading to 

Mu ^ qy 
Ml 



Ml), 



qui + 1 
Thus with M = {l + q)Mi, etc, (|) gives finally 

3 /. \ (3^-1) / ^ , X S+S/S/") 



P_ 

P 



l + q^ 

1 + g 



1 + qir] 



(8) 



(9) 



(10) 



where Pi is the initial period. 

As a check we note that if all the transferred mass is 
blown away from star 1, with the specific angular momentum 
of that star, we have i] — 0, f3 = 1; we take the limit of ( |l0[ ) 
as r; -> by noting that (1 -I- q-qf^'^ = exp[(3/»7) In(l + qrj)] 
and using I'Hopital's rule on the exponent to give the limit 
as 77 ^ as e^''. We then get 
3 /^ \ 2 



P_ 

P. 



1 + g, 
1 + g 



„3(q-q^) 



(11) 
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or P oc M2"^M"^exp[3M2/Mi], as found by e.g. King & 
Ritter (1999) for this case. 

We are interested in the minimum orbital period at- 
tained during the binary evolution. Regarding P in ( |lO| ) as 
a function of q, i.e. with Pi, qi fixed, we find 

1 dP _ 
P'd^ 



3/? - 1 3P + 5ri 



q 1 + q i + qv 

so that P is an extremum at 



(12) 



(l-r,) + v/(l-r?)^ + 9(r; + /j~^ 

' = '- = 3(. + /3-./3) ^''^ 

From ( |l2| ) we can show that 

1 d^'P 3 



+ 



(14) 



(3/3-l)(l-7?) 
P dq^ q^il + qv) ' (1 + q)^ {1 + qri) 

at q — qm, so P has a minimum there, assuming < j; < 1. 
Hence for given Pi,qi, the smallest orbital period which can 
be attained in evolution with mass loss is given by (|lO[) with 
q = q-m- For a given initial period, we still have the free- 
dom to vary the initial mass ratio qt. (This corresponds to 
the fact that the period Pi essentially determines the donor's 
mean density, almost independently of the primary mass. So 
a donor of given initial mass M2i may initiate mass trans- 
fer at various Pi, depending on its evolutionary state, with 
different qi, determined by M\i.) For fixed Pi, the minimum 
value of P clearly has a maximum, regarded as a function of 
qi, at qi = qm (see also Fig. H which shows the special case 
/3 = 1). In other words, A{P{q = gm))/(d(ji)[Pifixed] < 0. 
Since q decreases from its initial value qi through the evolu- 
tion this means that given an initial period Pi, the smallest 
value of the minimum period is given by the largest possible 
value qu oiqi. Hence given an initial period Pi, the minimum 
possible orbital period Pm{Pi) is given by 

qu \ l^ + qu \ / 1 + qmii 



(15) 
,9™/ \'i- + qmj yi + quvj 

with qm given by ([13|), and qu the largest possible value of 
the initial mass ratio qt. Although l3 is of course irrelevant 
in the conservative case {rj — 1), larger angular momentum 
losses in the wind allow even smaller ratios Pm/Pi for the 
same qi. 

It is likely that thermal-timescale mass transfer is still 
going on at the minimum orbital period, since by Roche 
geometry we have 

i?L oc /(g)M^/^P^/^ (16) 

where f{q) is a slowly increasing function of q given e.g. by 
the approximation from Eggleton (1983). Thus 



Rl _ g d/ A£_ 2P_ 
r£ ~ Jd^^ 3M '*' 3P' 



(17) 



which tells us immediately that the period reaches its miin- 
mum value at larger mass ratios than _Rl as g < 0, M < 0. 
In general -Rl initially shrinks much more rapidly than _Rto. 
Thus the condition R = _Rto signalling the end of thermal- 
timescale mass transfer is reached only once i?L is increas- 
ing after passing through its minimum value, i.e. thermal 
timescale mass transfer usually stops only a]ter the min- 
imum orbital period is reached. However it is conceivable 
that in some cases -Rtc shrinks very rapidly on mass loss, and 
individual cases must be checked. For the purposes of this 



paper we note than the true minimum period for thermal- 
timescale could if anything be longer than the values of Pm 
we find here. 

In what follows we will assume that /3 = 1, i.e. that the 
mass lost from the accretor has the same specific angular 
momentum as this star. This for example includes any form 
of mass loss from orbits with circular symmetry about the 
accretor, as might occur from an accretion disc. In this case 
([l5[) simplifies to 



p 



-qu 



-qm 



qmV 



quv 



5+3/v 



with 1 < Qm < 1.387 now given by 



qm = 



(l-'7) + V(l-^)' + 9 



(18) 



(19) 



and qu the largest possible value of the initial mass ratio qi. 
Figure ti shows the ratio Pm/Pi as a function of qi for the 
two extreme cases 77 = 1,0. 



3 MINIMUM PERIODS 

The work of the last Section, especially eqn (llSl), shows that 
to compute the ratio Pm/Pi of minimum to initial period 
we need to specify the largest possible value qu of the ini- 
tial mass ratio. For a white-dwarf accretor the short ob- 
served periods of RX J0537.7-7034 and IE 0035.4-7230 then 
constrain the mass-loss parameter rj. Clearly to reach such 
short periods it is preferable to start from the shortest ini- 
tial periods Pi, which in turn are given by assuming that the 
donor is still very close to the ZAMS when it initiates mass 
transfer. We can easily show that the initial donor mass M2i 
must be ^ 1.9Mq in order to explain the 3.5 hr period of 
RX J0537.7-7034, by iterating the formula (^, using the 
fact that we must clearly have qu > qm- For the lowest- 
mass white dwarf we consider {Mn = O.7M0) this implies 
M2i > O.7gmM0. Thus in the conservative case rj = 1 we 
have qm = l, M2, > O.7M0 , and Fig. | shows that Pj ^ 6 hr 
for a ZAMS star of this mass. We can now use Fig. H with the 
restriction Pm/Pi < 3.6/6 = 0.58 to find qu > 2.36 and thus 
M2i > I.66M0. Figure M now shows that Pi ~ 10 hr, and we 
may iterate using Fig. g to get q^ > 3.38, M2i > 2.37Mq. 
Further iteration fails, as Figured shows that the new higher 
estimate for M2i does not increase the estimate for Pi. In a 
similar way we find qu > 2.7, M2i > 1.89M0 for the other 
extreme case 77 = 0. 

This argument essentially fixes the minimum value of 
Pi at about 10 hr (see Figure ra), so the orbital shrinkage 
required to explain the 3.5 hr period of RX J0537. 7-7034 
has Pm/Pi < 0.35. 

Given an initial white dwarf mass Mu and mass ratio 
qu we can predict the period ratio Pm/Pi for any given rj. 
The initial mass ratio qu must obey two constraints (added 
as vertical lines in Fig. B): 

(i) the system must avoid the 'delayed dynamical in- 
stability' (DDI; see Webbink, 1977; Hjellming, 1989), which 
occurs when sustained thermal-timescale mass transfer ex- 
poses inner layers with a flat entropy gradient, and 

(ii) the white dwarf mass cannot exceed Mc before the 
minimum period Pm is reached. 

In practice the first of these constraints requires qu < 
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3 

qi 



Figure 2. Pm/Pi as a function of qi for r; = 1, according to (|lq). 
Additional constraints for various initial WD masses resulting 
from the Chandrasekhar limit are indicated by dashed lines in 
the conservative case. They mark the critical values of qi where 
the WD mass has grown to Mc precisely upon reaching Pm- An 
upper limit at g^ = 3.0 due to delayed dynamical instability (DDI) 
is also indicated. 
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Figure 3. Typical range of Pi versus M2i for different WD 
masses, where the lines indicate the main sequence boundaries. 
For initial donor star masses in the expected range between 
1 . . . 5Mq the short period edge (= ZAMS) lies at about 10 hr, 
and is fairly constant for M2i ft 1.3Mq . The additional line marks 
the terminal main sequence (TMS), more precisely the point of 
minimum T^g (where applicable). 



gDDi — 3. Hjellming (1989) found this value for a donor 
with M2i = SMq near the terminal main sequence (TMS), 
whereas Kalogera & Webbink (1996) seem to prefer even 
smaller limiting ratios around 2.5. Kolb et al. (2000) used 
Mazzitelli's stellar code to calculate a test sequence with 
constant primary mass O.75M0 where mass transfer starts 
from a SMq near-TMS star. A DDI occured at donor mass 
2.6M0, in perfect agreement with Hjellming's prediction. 
The critical maximum mass ratio goDi depends on the stel- 
lar structure, and therefore on the stellar input physics. In 
particular, ^ddi is probably sensitive to the degree of con- 
vective core overshooting, as this determines the size of the 
convective core. This is highlighted by the fact that Kolb 
et al. (2000) find qnui = 2.9 for early massive case B mass 
transfer, while Tauris et al. (2000) find qooi = 3.6, using an 
updated Eggleton code. 

The second constraint is quite severe for large initial 
white dwarf masses and for 77 ~ 1. In particular, for Adu = 



0.7Mq it requires qn < 3.11, in mild contradiction with the 
requirement qn > 3.38 we found above. 

Figure W shows Pm/Pi versus r) for various values of qn 
The stars denote combinations Mii,qii,rj where Mi reaches 
Mc precisely at q = qm- To the right of these positions we 
take Pm as the minimum period actually achieved, i.e. the 
period where the white dwarf reaches Mc- 

Each panel of this Figure shows: 

1. The horizontal line Pm/Pi ~ 0.35, i.e. the upper limit 
required by RX J0537.7-7034. 

2. Only those curves which actually manage to cross this 
line for 77 > 0. In particular we plot the curve for the limiting 
value of qi such that the curve just crosses Pm/Pi = 0.35 at 
rj — 0. This gives a lower limit on qi which is the same for 
each mass of the WD. 

3. At some of the extreme solutions, i.e. wherever the 
line Pm/Pi = 0.35 is crossed before a SN occurs, the current 
values of Mi , M2 are given. Note that with growing Adu the 
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Figure 4. Pm/Pi versus rj for various g^ as labelled. Only the 
part of each curve in thick linestyle reaches Pm before the WD 
grows to the Chandrasekhar mass and a SN la occurs, with an 
asterisk marking the critical 77. Beyond that the current period 
at the SN event has been given as the smallest possible fraction 
Pm/Pi during the thermal— timescale mass transfer phase. The 
horizontal dotted line at Pm/Pi = 0.35 marks the fraction at least 
required to create a system like RX J0537.7— 7034 starting from 
the ZAMS with Pi ^ 10 hr. The numbers in brackets show pairs of 
masses (Mi((jm), M2{qm)) each belonging to a set of parameters 
(77, qi) just fulfilling this requirement and marked by a diamond. 
DDI constraints (that would exclude the two larger initial mass 
ratios) have not been considered in this graph. 



donor mass has to go up as well for the same qi. Furthermore 
Ml at the minimum period comes closer and closer to the 
Chandrasekhar mass as we increase r). 



4 DISCUSSION 

Figure shows already that the short orbital period of 
RX J0537.7-7034 poses very severe constraints if this system 
results from thermal-timescale mass transfer. In particular 
(a) Conservative evolution {jj — 1) is possible only for 
initial mass ratios qi > 3 which probably make the system 
vulnerable to the delayed dynamical instability. Full evolu- 
tionary calculations are required to check if there are any 
evolutionary tracks which avoid it. If such tracks exist, the 
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Figure 5. Borderline in the plane of r) and qi needed to obtain 
Pm/Pi = 0.35. The allowed range is limited to the small area 
above the curved line and below the horizontal DDI limit at qi = 
3.0. Additional constraints are shown for various Mu (dashed 
lines). The upper limits to gi(r;) (or »y(gi)) result from the WD 
growing beyond its Chandrasekhar mass before reaching Pm, see 
Fig. W and text for details. Initially massive WDs with ^ I.QMq 
restrict the allowed range to 77 ^ 0.3 independently of DDI. 



initial white dwarf mass in RX J0537. 7-7034 must have been 
low {^ O.7M0), but the current system masses must be fairly 



high (e.g. qi = 4 requires Mi = 1.17Mq,M2 



2.33Mc 



These are of course in conflict with the mass estimates of 
Greiner et al. (2000), but this may not by itself be fatal (see 
the remarks in the Introduction). 

(b) Non-conservative evolution (77 < 1) does allow 
tracks with g; < 3 which probably avoid the delayed dy- 
namical instability. However the ranges of 77 and qi are still 
very tightly constrained. Figure a shows that 2.7 f~, qi ^ 3, 
< ?7 ;$ 0.3, with 77 and qi correlated as in Figure H. 
For Mu ^ I.OM0 77 has to be even lower, and generally 
upper limits on qi are given for each Mu in addition to 
DDI. Note that for the largest values of 77 the white dwarf 
in RX J0537. 7-7034 is predicted to be close to the Chan- 
drasekhar limit unless the DDI limit imposes a more severe 
restriction than the lowest possible WD mass (as is the case 
in Fig. I for g, ^ 3, Mu ^ 0.7Mq). 

Even given these tight constraints, the evolutions dis- 
cussed above require the system to have come into con- 
tact with the donor still very close to the ZAMS. As or- 
bital angular momentum losses (e.g. via magnetic brak- 
ing) are probably negligible for the likely initial masses 
M2i ^ 2Mq, this also requires the initial orbital separa- 
tion to lie in an extremely narrow range. Evidently to make 
a system like RX J0537.7-7034 or IE 0035.4-7230 by the 
thermal-timescale route is a very rare event. In line with 
this, Deutschmann (1998) found no orbital periods shorter 
than about 6 hr in his detailed calculations with solar metal- 



© 0000 RAS, MNRAS 000, 000-000 



6 A.R. King et al. 



licity. Furthermore these requirements have been derived as- 
suming that RX J0537.7-7034 and IE 0035.4-7230 are just 
at the minimum period Pm, so the conditions might be even 
harder to meet. The probability of observing such a system is 
larger near the end of the thermally unstable phase because 
the mass transfer rate decreases (in both Deutschmann's and 
our own full computations). 

We have so far neglected the effect of tidal interactions 
on the orbital evolution of the binary. Tauris & Savonije 
(2000) show that these can be important in low-mass X- 
ray binaries by translating spin angular momentum losses 
(via e.g. magnetic stellar wind braking) into orbital losses 
when tidal synchronization occurs. However for the binaries 
we consider here, any effect before the beginning of mass 
transfer simply allows a shorter Pi for a wider range of sys- 
tems. Once synchronism is achieved, the angular momentum 
of a lobe-filling donor is less than about 10 per cent of the 
orbital angular momentum (because the gyration radii of 
typical donor stars is small, r"^ < 0.10. . .0.20). Even trans- 
ferring all of this to the orbit in the course of the evolution 
would lead to only a marginal shift towards longer minimum 
periods, leaving our conclusions unchanged. 

The work of this paper suggests that while the thermal- 
timescale mass transfer model for the supersoft X-ray bina- 
ries hcis many desirable features, it may not be possible to 
use it to describe all of the supersoft binaries, as well as 
SNe la progenitors. For example, if highly non-conservative 
mass transfer is as common as seems to be required to ex- 
plain RX J0537.7-7034 or IE 0035.4-7230, this would make 
building up the white dwarf mass to produce a Type la su- 
pernova highly problematical. It may therefore be necessary 
to consider other possibilities (e.g. van Teeseling & King, 
1998; King & van Teeseling, 1998). 
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